Polyelectrolyte Layer-by-Layer Assembly from Self-Consistent Field Calculations 
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We have modeled the layer-by-layer assembly process of flexible polyelectrolytes on flat surfaces. 
The multilayer has a three-zone structure. An exponential growth is found for the first several layers, 
followed by a linear growth for subsequent layers evolving towards a steady state. While adjacent 
layers are highly interpenetrating, stratification can be seen for every four or more layers. The 
effects of surface charge density, bulk salt concentration, and solvent quality on the thickness and 
internal structure of the multilayer are also studied. Our results agree with experimental findings. 

PACS numbers: 82.35.Rs, 68.65.Ac, 68.55.Jk 
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Proposed in early 90 's, the polyelectrolyte (PE) layer- 
by-layer assembly technique has attracted exponentially 
increasing interest due to its simplicity, versatility, and 
potential applications. 0,0,13 Through alternating expo- 
sure of a charged substrate to solutions of polyanions and 
polycations, hundreds of thin adsorbed layers of polyions 
can be readily built-up on the substrate. In addition to 
charged polymers, charged nano-particles and platelets 
have also been used.Q This technique has a wide variety 
of potential applications, including surface modification, 
sensors, conducting or light-emitting devices, drug deliv- 
ery, nano-reactors, etc.0 

Our understanding on the formation mechanism, inter- 
nal structure, and molecular property of PE multilayer, 
however, is still at an early stage. Q In great contrast to 
thousands of experimental papers on the layer-by-layer 
assembly, very few theoretical 0,0 or simulation^ 0,0 
studies have been reported. In general, the build-up 
of PE multilayers is based on the charge inversion, i.e., 
the charge of newly adsorbed PE inverts (overcompen- 
sates) the existing charge in the him (including the bare 
substrate charge). Using a self-consistent field (SCF) 
theory 0, we recently examined in detail the effects of 
various parameters in the adsorption of flexible PE on 
flat surfaces; strong charge inversion was found for rela- 
tively long PE on oppositely charged, attractive surfaces 
in poor solvent at high salt concentrations. jlO| In this 
work we use the same formalism to model the process of 
PE layer-by-layer assembly and to examine the thickness 
and internal structure of the multilayer. 

We only consider monovalent systems, and assume that 
all ions carrying the same type of charge are identical. We 
further ignore the volume of small ions, and impose the 
incompressibility constraint everywhere in the system. 
Under the ground-state dominance approximation^^! 
the SCF equations describing the adsorption of charged 
polymer A, in the presence of polymers 1 and 2 (of which 
density profiles are fixed) , from its bulk solution of small- 
molecule solvent S (having the same density po as all the 



polymer segments) onto a flat surface are given bv|lfjj 
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where x denotes the direction perpendicular to the sur- 
face (placed at x — 0); <fip{x) is the normalized (by po) 
segmental density of polymer P (= A, 1, and 2); ip(x) 
is the electrostatic potential (in units of fc^T/e, where 
kg is the Boltzmann constant, T the absolute tempera- 
ture, and e the elementary charge), and we set UbT = 1 
hereafter; uj a (x) = N{J2 p (xpa - xps)0p(x) + Xas[1 - 

£p - M 1 ~ Ep^pWH. where XPP' ( p ' = A > 

1, and 2) and xps are the Flory-Huggins parameters 
(xpp = 0); gp(x) = —vpppip(x) for the smeared charge 
distribution where each segment of polymer P carries 
charge vpppe (a model for strongly dissociating PE), and 
gp(x) = ln[l — pp + pp exp(— vpip(x))} for the annealed 
charge distribution where each segment of P has a prob- 
ability pp of carrying charge vpe (a model for weakly 
dissociating PE), and we set V\ = —v% = — 1; 0a, b an d 
c Sj b are the normalized (by po) segmental density of A 
and salt concentration (we assume that one salt molecule 
generates one cation), respectively, in the bulk solution 
placed at x = I > 0. N ote t hat we have normalized 
the distance by R g = a \J N/6, where a is the statisti- 
cal segment length of polymers (assumed to be the same 
for all P), and is an (arbitrary) chain length chosen 
for the normalization; we have also normalized the di- 
electric constant e (assumed to be uniform for x > 0) 
by 27r poe 2 a 2 / (3k bT), where e before the normalization 
is in units of 47re with e = 8.854 x 10" 12 (A-s) 2 /(J-m) 
being the permittivity of vacuum. Note further that we 
have chosen = and </>a(0 = 4>A.b, which effectively 
couple the bulk solution to our system for both ions and 
polymer A, and thus give ujq — tu A (l); they are satisfied 
by choosing large enough system size I such that both 
\4>a(1) — <t>A,b\ and |^>(i)| are less than 10~ 10 . 

The boundary conditions for the mean-field equations 
are lOj: (din \/4>\/dx)\x=o — d -1 , where d" 1 repre- 
sents the short-range (non-Coulombic) surface-polymer 
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FIG. 1: (Color online) Polymer segmental density profile of 
the i th layer 4'a > w( x )- The layer number i is marked near the 
corresponding profile. <tsf = 0.1, c S) j, = 0.05, and X12 = 0. 

interactions and includes both energetic and entropic 
contributions; {d-\/(f>A/dx)\ x —i — 0; (dip/dx)\ x =o = 
— asFV^V /e, where the normalized surface charge density 
<jsf = PaVS/ipo®) > 0, and Co is the number of surface 
charges (in units of e) per unit area; and (dip / dx)\ x= i = 
0, which ensures the overall charge neutrality of our sys- 
tem (including the surface charge) . The mean- field equa- 
tions are solved by a relaxation method (l2^. 

We do the following to mimic the process of layer-by- 
layer assembly of polymers 1 (polyanions) and 2 (poly- 
cations): (1) For the i th deposition, A is set to represent 
polymer 1 when i is odd and 2 when i is even. After 
solving the mean-field equations, the amount of adsorbed 
polymers in this deposition, rW = y/~N Jq[4>a( x ) ~ 
4>A,b]dx, is calculated. (2) A washing step follows each 
deposition and modifies the polymer density profile from 
4>^{ x ) to <f$ w {x). Since it is not possible to distin- 
guish between adsorbed and unabsorbed chains under the 
ground-state dominance approximation, we (somewhat 

arbitrarily) choose 4$ w (x) = $ (x) exp[tx 2 /(x — Xtu )] 

for x < Xw and 4>^ w (x) = for x > Xw\ where the 

determination of the "washing point" Xw G (0, 1] is ex- 
plained in Ref. [l^, and t > is calculated according 

to VN ff™ <f>^ >w (x)dx = rW. (3) 0i (x) and fa(x) to be 
used in the i + 1 deposition are accumulated for polyan- 
ions and polycations, respectively, from 4>x w {x) for all 
j < i. That in each deposition <j)\(x) and ^(x) are 
fixed represents the irreversibility of PE adsorption Q; 
we therefore model the multilayer growth using a series 
of kinetically trapped states. 

In this Letter we study the multilayer formation on 
an indifferent surface (d _1 = 0) using symmetric PE 
(pp = 0.5 and 0a, b — 7.5 x 10~ 4 ) both having smeared 
charge distribution in a poor solvent (xps = 
Fig. n shows the polymer segmental density profile of 




FIG. 2: (Color online) Total segmental density profiles of 
polyanions and polycations after washing, 4>i{x) and <f>2 (x), 
respectively, and the electrostatic potential profile before 
washing, ip(x), after the 60 th deposition. The inset shows the 
same profiles after the 6 th deposition, ctsf = 0.1, c s fi = 0.05, 
and X12 = 0. 



the i th deposition after washing (referred to as the i th 
layer), 4*^ w {x), on a positively charged substrate of 
(7 g p = 0.1, where a bulk salt concentration of c St b — 0.05 
is used in each deposition. [TH Under such conditions, 
the first layer strongly overcompensates the substrate 
charge: Ect^/^sf ~ —1-36, where Eer^ = o^f + 
J2]=i VjppjirV) is the overall charge of the multilayer 
(including the substrate charge) after the z th deposition. 
We set X12 = in this case; it is therefore the presence 
of the first layer of polyanions that leads to the adsorp- 
tion of the second layer (i.e., T^ 2 ' > 0), which carries 
the same type of charge as the substrate. Similarly, we 
have rw > r('~ 2 ) for i = 3 to 8. Subsequent layers 

(81 

have similar density profiles to (f> A w (x), and are further 
away from the substrate; consistent with experimental 
observations El- Fig. ^ also shows that, while adja- 
cent layers are highly interpenetrating, stratification can 
be seen for every four or more layers. 

The total segmental density profiles of polyanions and 
polycations, <j)\{x) and 0a {x), respectively, after the 60 th 
deposition and washing are shown in Fig. |2 also shown 
is the electrostatic potential profile tp(x) after the 60 th 
deposition (before washing). Three zones can be identi- 
fied in the multilayer: In the interior of the film (Zone 
II), 4>\{x) and 4>2{x) are almost constant and ip(x) w 0. 
Here, <pi(x) ~ 4>2{x) due to the symmetric PE we used, 
and the concentrations of ions (before washing) are close 
to those in the bulk solution. We see small oscillations 
in 4>\(x) and 4>i{x) (which are out-of-phase), indicating 
the weakly spatial enrichment of polyanions and polyca- 
tions; ip(x) also oscillates, in-phase with ^>s(x), in Zone 
II. Close to the substrate (Zone I), <f>i{x) > <f>2(x) due to 
the electrostatic interactions between the substrate and 
polymers. In the immediate vicinity of the substrate, 
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TABLE I: Average total polymer segmental density in Zone 
II of the multilayer 0pbm- Since the total polymer segmental 
density slightly varies in Zone II, we use three times its stan- 
dard deviation as the error, ctsf = 0.1 and c Sl (, = 0.05. Refer 
to text for more details. 
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Layer i 

FIG. 3: (Color online) Semi-logarithmic plot of the overall 
charge of the multilayer Ect^' as a function of layer number i 
(Ect w < when i is odd and Ect w > when i is even). The 
straight lines are fits for the first several layers exhibiting an 
exponential growth. X12 = 0. 

ip{x) is positive and larger than that in Zone II. At the 
outer boundary of the multilayer (Zone III), both <fri(x) 
and <p2 (%) decay to 0, and ip(x) exhibits a large peak (pos- 
itive when i is even and negative when i is odd). Such 
a three-zone structure was first proposed by Castelnovo 
and Joanny for multilayers of flexible and symmetric PE 
on an indifferent surface in 0-solvent at high salt con- 
centrations, where the polymer density in Zone II was 
assumed to be constant and equal to that of neutral 
complex formed in the bulk.p| It is also supported by 
the molecular dynamics simulations of Patel et al, who 
studied the process of layer-by-layer assembly of symmet- 
ric PE on a planar substrate in a poor solvent with no 
added salt; the density oscillations in Zone II were also 
observed by these authors. Q 

For a homogeneous system, the SCF theory reduces 
to the Flory-Huggins theory modified for PE[T(j; in our 
case of symmetric PE 1 and 2 (with chain lengths Np — > 
oo ) having smeared charge distribution in solvent S, the 
free energy (of mixing) per polymer segment or solvent 
molecule becomes 

f=^hx^ + (l-<P) ln(l - 4>) + Xps^(1 -<t>) + Xi2^ 
+ (pp0 + 2p.)ta**±^- (3) 
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Layer i 

FIG. 4: Overall thickness of the multilayer T,D^ as a func- 
tion of layer number i. The inset shows the thickness of the 
60 th layer Z)' 60 ' as a function of bulk salt concentration c a ^- 
X12 = 0. 

where 4> (— 2</>i = 2<f>2) and c a are the total polymer 
segmental density and salt concentration in the homoge- 
neous system, respectively. From this, the binodal curve 
describing the phase equilibrium between a polymer-rich 
phase (where the total polymer segmental density is de- 
noted by C mp) and a polymer-poor phase (where the 
polymer density is 0) can be derived, which is given by 
ln(l - 0cm P )/0cmp + (xps - Xi2/4)0cm P + 1 = 0. We also 
define m in as the total polymer segmental density that 
minimizes / (with c s = c S] b), and <?!> m in.N as the one that 
minimizes the neutral part of /, i.e., the first line of Eq.0 
Table H] shows that the average total polymer segmental 
density in Zone II of the multilayer, denoted by </>pemj 
is very close to </> m in.Ni instead of (/> cmp ; our results also 
show that 4>pem is independent of c s ^ and ctsf- 

The inset of Fig.|21shows that Zone I is almost fully de- 
veloped after the first six layers; subsequent layers start 
to move away from the substrate (as shown in Fig. 
The first six layers exhibit an (approximately) exponen- 
tial growth, different from the subsequent layers evolv- 
ing towards a steady state in a damped oscillatory man- 
ner. This is shown in Fig. [3J where EcrW < when i 
is odd and EcrW > when i is even; each layer there- 
fore inverts the overall charge of the multilayer. An 
almost constant |Ecr^^| is reached in the steady state, 
which implies that the charge of the newly adsorbed 
layer i fully inverts the overall charge of the existing 
multilayer, i.e., Vj^paT® fEa^ 1 ) « —2; this was also 
found by Dobrynin and co-workers 0, @- Fig. 0] shows 
how the overall thickness of the multilayer, £Z)M = 

2^j2*E5=i^A,U s ) diB /£5=i r0 ' ) ' evolves with the 

layer number i; an (approximately) linear growth is seen 
after the first six layers. These are in good agreement 
with experimental observations^,!^. 

The above characteristics are seen in all of our calcula- 
tions when multilayer forms. The effects of ctsf and c S) b 
are also shown in Figs. |21 and 0] The surface charge den- 
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sity (tsf affects only the exponential growth of the multi- 
layer (and thus ££)W), but not the adsorbed amount J?w 
(or equivalently EerM) and layer thickness in the steady 
state. The bulk salt concentration c s ^ plays a critical 
role in determining and the layer thickness (and thus 
£Z)W); the inset of Fig.0]shows a quadratic relationship 
between the thickness of the 60 th layer (where a linear 
growth is reached in all cases), L>( 60 ) = EL>( 60 ) - SZ)( 59 ), 
and c Sj b- Note that our mean-field theory does not give 
strong charge inversion an d multilayer formation at 
low salt concentrations, since the fluctuations and lateral 
correlations in the system are ignored. 

Finally, we also used 0-solvent where xps = 0.5 (with 
usf = 0.1 and X12 = 0); no multilayer formation is seen, 
even at c s .b = 0.1 where the charge inversion of the first 



layer (EcrW / <tsf ~ —0.46) is stronger than that in the 
poor solvent at c s ,b — 0.01 (where Ect^/csf ~ —0.27 
and multilayer forms). As noted above, > r^ -2 ^ is 
observed during the exponential growth of the multilayer 
in the poor solvent (except for i = 3 in the c s .b = 0.01 
case). The opposite, however, occurs in the ^-solvent, 
where monotonically decays towards 0. The solvent 
quality therefore also plays an important role in the mul- 
tilayer formation. This agrees with the molecular dynam- 
ics simulations of Dobrynin and co-workers, who found 
that multilayer forms under poor solvent conditions but 
not under good solvent conditions 0, 0- 
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